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In this paper, we do three things in the study of deuteron break-up by high energy neutrino beams. (1) 
We present previously unpublished data on neutrino induced backward protons from deuteron targets; (2) 
we calculate the contributions from both the two-nucleon (2N) and six-quark (6q) deuteron components, 
which depend upon the overall normalization of the part that is 6q; and (3) we suggest other signatures 
for distinguishing the 2N and 6q clusters. We conclude that the 6q cluster easily explains the shape of the 
high momentum backward proton spectrum, and its size is nicely explained if the amount of 6q is one or 
a few percent by normalization of the deuteron. There is a crossover, above which the 6q contribution is 
important or dominant, at 300-400 MeV/c backward proton momentum. 
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I. INTRODUCTION 

It is worthwhile to examine and probe the deuteron 
as a target from many viewpoints because of the role 
it plays as our main source of information on the neu- 
tron. This has become of increasing interest in recent 
years with unexpected results being found for the nu- 
merous sum rules in which neutron structure functions 
enter. These include the results for the Ellis- Jaffe sum 
rules |jl|,||, for the Gottfried sum rule |^, and for the 
Bjorken sum rule MM, prompting us to re-examine as- 
sumptions made in most analyses of deuteron data to 
extract neutron information. The deuteron is more than 
the sum of its proton and neutron parts and since experi- 
ments on free neutrons can not be done, this nonadditive 
portion is important to evaluate by any means available. 
We believe, in fact, that earlier experiments [bIq] make it 
clear that the deuteron is the simplest nucleus exhibiting 
the EMC effect which, in turn, must affect the extraction 
of the neutron structure function that enters in these sum 
rules [^. 

In the theoretical part of the present paper, we enlarge 
upon the work we introduced earlier M and apply it to 
some old data and to new, previously unpublished data. 
We have studied deep inelastic v and D scattering from 
various targets, focusing on reactions that produce high 
momentum backward protons. Backward means relative 
to the incoming neutrino or antineutrino and high mo- 
mentum means relative to kinematic limits upon back- 
ward momentum imposed in terms of quantities such as 
the target mass and momentum fraction carried by the 
struck quark. 

We concluded, based on deuteron target data ll^l ob- 
tained at CERN with v and D beams, that the momen- 



tum distribution of the backward protons was consistent 
with production from a multiquark component (6q) in 
the deuteron and was difficult to explain if produced 
via break-up of a two nuclcon (2N) quantum mechani- 
cal state with a simple conventional wave function. The 
particular problem with the latter picture, where the neu- 
tron and proton substantially maintain their character as 
nucleons pT| ] , was that there were more large momentum 
backward protons than expected from typical neutron- 
proton wave functions. However, in our previous work 
we were only able to calculate the shape of the backward 
proton spectrum in the two cases and not the absolute 
normalization or even the relative normalization of the 
6q and 2N contributions. 

In this paper, we will present normalized calculations 
for backward proton production in deep inelastic exper- 
iments for both the 2N and 6q deuteron components. 
The 6q calculation includes factors of the fraction of the 
deuteron that is 6q and of the fragmentation rate of the 
residue of the 6q cluster into protons. A numerical un- 
certainty in the latter rate for high momentum protons 
leads to a factor circa 2 uncertainty in fixing the fraction 
of 6q state in the deuteron; otherwise the calculation is 
well determined. 

We denote the probability of finding the 6q configura- 
tion in the deuteron as /, with / expected to be between 
0.01 and 0.07. Possible values for / have been calculated 
from deuteron wave functions as the probability for the 
nucleons to overlap, and Sato et al. ||l2| claim a value for 
/ in the middle of this range for a typical nucleon-nucleon 
potential. The largest deuteron probability, / = 0.074, 
was used [Q in describing high energy SLAC electron- 
deuteron data at x > 1.0. Also, studies of the deuteron 
electromagnetic structure functions have been used to es- 



timate / to be a few percent [|Ij] . We here will find that 
the 6q contribution in deep inelastic scattering, even with 
/ only one or two percent, can be quite large for energetic 
backward protons. Throughout, although we are mind- 
ful of many possibilities (see e.g., Ilq|), we simplify by 
speaking of the short range baryon number correlation 
as either 6q or 2N. Within this limitation, we can further 
say that for most of the conventional 2N wave functions, 
the 6q state not only can but also must contribute a ma- 
jor share of the cross section for energetic backward pro- 
tons. We will show that at 500 MeV/c backward proton 
momentum the modern nucleon-nucleon potential that 
comes closest to our data needs at least a 60% additional 
contribution. 

On the experimental side, we shall in the next Sec- 
tion discuss a Fermilab neutrino-deuterium 15-ft bub- 
ble chamber experiment (E545) which has obtained data 
(previously unpublished) on backward production of pro- 
tons. Also, we shall comment on some published exper- 
iments which likewise measured high-momentum back- 
ward proton events. 

In Sect. Ill, we give the theoretical formulae used in 



calculation of backward proton spectra to compare with 
the Fermilab and CERN deuteron break-up cross sections 
for backward protons. A comparison of the new data 
with previously published Argonne, Brookhaven, CERN, 
and Fermilab backward proton production data is also 
given. The comparisons of the prediction given by mod- 



IV. The final 



ern potentials with our data is in Sect. 

Section is devoted to further discussion and presentation 

of conclusions. 



II. FERMILAB E545 NEUTRINO-DEUTERIUM 
DATA 

Previous analyses of the E545 neutrino-deuterium data 
have concentrated on either extracting vn interaction 
data, or on the distributions of lepton or hadron vari- 
ables in vd scattering. A typical analysis would sepa- 
rate the data into even-prong 'Vn" and odd-prong 'Vp" 
events, where a visible proton spectator is ignored in the 
prong count. The observed 'Vn" events were assumed 
to be a sample of vn interactions depleted by "rescatter- 
ing" within the deuteron nucleus. The "rescattered" lyn 
events, in turn, would appear in the 'Vp" event sample. 
In extraction of ratios of vn to vp cross sections the frac- 
tion of such rescattered events was estimated to be in the 
6-12% range. 

An excess of high momentum proton spectators, com- 
pared to standard deuteron wave function predictions, 
is known to be associated with the "i/n" event sample. 
When commented on in previous analyses, it would gen- 
erally be noted that the origin of this excess is unknown, 
but given that the excess accounts for less than 1% of 
the 'Vn" events any effect on the overall distributions or 
measured quantities would be negligible. In the present 



analysis, we explicitly examine the vd target fragments 
in the E545 data, and compare the distribution of back- 
ward protons with a deuteron wave function plus a small 
six-quark component. 

Many details of the E545 experiment have been pub- 
lished ijl^l. The data discussed here were previously 
presented at an APS meeting [|^, where the empha- 
sis was on establishing the existence of a "rescattering" 
phenomena which depleted the observed spectator pro- 
ton (neutron target) event sample, and in estimating its 
frequency. 

The E545 data are from a 320,000 frame exposure of 
the deuterium- filled Fermilab 15-ft bubble chamber to a 
wide-band single-horn focused neutrino beam produced 
by 4.8 X 10"'^^ 350 GeV/c protons incident on a beryllium 
oxide target. The anti- neutrino component of the beam is 
« 14%. The film was scanned twice, and events with two 
or more charged tracks produced by incident neutral par- 
ticles in a 15.6 m'^ fiducial volume were accepted for anal- 
ysis. All charged tracks were digitized and geometrically 
reconstructed. Topology-dependent weights are applied 
to the data to compensate for scanning and processing 
losses and for those events failing geometric reconstruc- 
tion. The average processing times scanning efficiency is 
0.80. Cuts are applied to the two-prong events to remove 
K'^ and A decays and 7 conversions from the data. 

A kinematic technique which uses only the measured 
momenta of the charged particles is used to select a 
sample of charge current events. Only events for which 
J2pl > 5 GeV/c, where pl is the component of labo- 
ratory momentum in the beam direction and the sum is 
taken over all charged particles, are included in the anal- 
ysis. The muon candidate is identified as that negative 
track in the event with the largest component of momen- 
tum transverse to the incident neutrino direction. Those 
events for which the component of the /i^ candidate's 
momentum transverse to the vector sum of the momenta 
of the other charged particles in the event is greater than 
1.0 GeV/c are accepted as charge current events. 

The incident neutrino energy of the selected charge 
current events is estimated using transverse momentum 
balance: Ei, — p'f + p^ -\- |py -I- pf}\p^ /pi}, where the 
symbols p^ and p refer to the muon momentum and the 
vector sum of the charged hadron momenta, respectively. 
Only events with Ey > 10 GeV are accepted for analysis. 

A Monte Carlo simulation indicates that the sample se- 
lected according to the above criteria includes 79% of the 
vd charge current events, with the fi~ correctly identified 
in 98% of the cases, and with a 3% contamination due to 
vd neutral-current events and 1% due to Pd events. 

The corrected number of i>d events in the sample is 
15,129, with an average neutrino energy (Ei,) — 50 GeV. 
Of these events, 459 have an identified proton with mo- 
mentum magnitude greater than 160 MeV/c whose di- 
rection is backward with respect to the incident neutrino 
direction. (Significant visibility losses occur for protons 
with momentum less than 160 MeV/c, and hence are 
not presented.) The identity of the backward protons 



was verified by re-examining all such tracks on the scan 
table. The momentum distribution of the backward pro- 
tons is given in Table |. This data will be discussed in 
Sect. 3, together with the proton spectrum from a vd 
and Dd exposure of BEBC by the WA25 collaboration at 
CERN 0. 

The published Fermilab E545 spectator proton spec- 
tra from quasi-elastic I'd scattering {vd — > ji^pps) |lq ] 
will also be discussed, together with similar distributions 
from vd bubble chamber experiments at Brookhaven (80- 
in) and Argonne (12-ft) |§. 



III. THEORETICAL DISCUSSION 

We will give the expressions for the charge current in- 
clusive cross sections of neutrinos hitting a deuteron and 
producing a backward proton, ps, a forward lepton, ^~, 
and anything else, AT, 
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e +PB + X, 



(1) 



in both the 2N and 6q models. We will see that the 
shape of the backward proton spectrum is different in the 
two cases, and will see if the calculations can match the 
data. The two contributions arc not mutually exclusive, 
and we shall add them incoherently, weighting the 6q 
contribution by fraction / and the 2N contribution by 
(1 — /). We expect that / will be on the order of a few 
percent, but that none-the-less the 6q contribution could 
be large at large backward proton momenta. 

The fully differential cross section is differential in x, y, 
a, and px, which are the experimentally measurable vari- 
ables. These variables are the struck quark momentum 
fraction. 



the fractional lepton energy loss, 
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the light front momentum fraction of the backward pro- 
ton, 
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(with pz defined positive for backward protons); and the 
transverse momentum of the proton relative to the direc- 
tion of the incident neutrino, pt- 

For the 2N model, using quark distribution functions 
appropriate to describe striking the neutron, with the 
argument changed from a: to ^ = 2:/ (2 — a) because the 
neutron is moving, and having the proton emerge with 
probability given in terms of the deuteron wave function, 
we find the cross section 
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where _D„ (^) is ^ times the distribution function of down 
quarks in the neutron (n), etc., 7 = Ep/mN, Oq is the 
point Fermi weak interaction (with strenght Gp) cross 
section, 
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and i/j is the wave function of the deuteron normalized 
by 



dad^PT \4'{a,PT)\ 
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The corresponding cross section for the 6q component 
of the target can be written in terms of the probability 
distribution of a quark in the 6q cluster and in terms of 
the probability, D^/^q , for the residuum of the 6q state to 
fragment into the proton. This time, since the deuteron 
or 6q cluster is stationary in the lab, x is directly — in 
the scaling limit — the momentum fraction of the struck 
quark. We have 
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:Dp/5q{z,px), (8) 



where we have included just Dq, the down quark distri- 
bution for the 6q cluster times x, on the grounds that 
we will need the literal 5q residuum (which comes from 
the 6q Fock component of the nominal 6q cluster) to get 
the highest momenta backward protons. The first argu- 
ment of the fragmentation function is the light front mo- 
mentum fraction of the proton relative to the five quark 
residuum, or 



Momentum range (MeV/c) 


Number of events 


160-200 


187 


200-240 


100 


240-280 


61 


280-320 


37 


320-360 


30 


360-400 


14 


400-440 


10 


440-480 


10 


480-520 


9 


520-560 


1 


560-600 






TABLE I. Momentum distribution of backward protons in 
15,129 vd charge current events from Fermilab experiment 
E545. 



Presently reported data on the backward proton mo- 
mentum spectrum uses protons gathered from the entire 
backward hemisphere. Hence, we too will integrate over 
the backward hemisphere, to allow direct comparison to 
experiment. We will integrate over x and y also. The 
term with explicit y dependence gives a small contribu- 
tion. Then, 
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where p = \p\ is the backward proton momentum (and 
we used ada/da d^pr — Edcr/d^p) . What we want is the 
weighted sum of the 2N and 6q contributions. 
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The plotted curves are based on the above formulas 
plus some choices for the deuteron wave function, the 
quark distribution functions, and the fragmentation func- 
tion of the 5q residuum. 

The wave function will be a light front wave function. 
(A simple use of a non-relativistic wave function conflicts 
with the kinematic bound on maximum backward proton 
momentum.) It is related to non-relativistic wave func- 
tions by 



|V'(a,PT)l^ = \i}LF{a,PT)\ 
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where the arguments of the non-relativistic wave function 
are obtained from 
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The normalization is 



d^k \ipNR{kz,kT)\ = 1 



and the factor above comes from the Jacobian in 
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We use several different deuteron wave functions, but 
start with a Hulthen wave function, which is still in com- 
mon use |^,|l|,|l|. 
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The quark distributions for the nucleon are the set 
CTEQIL p3]. (Some old and simple quark distribu- 
tions HI] give results about the same.) For the 6q clus- 
ter, we use the Lassila-Sukhatme model "B" quark dis- 
tributions |24|. These distributions are based on quark 
counting rules and physical logic and describe the EMC 
data. Models "A" and "C" are not very different for the 
present purposes and are omitted from the figures mainly 
to avoid clutter. Model "B" has 



D,{x) = 3xl.85J^{l-^y\ 



(20) 



The fragmentation function for the 5q residuum is 
taken in a factorized form, 

(21) 

where A — 310 MeV. The spectrum of protons for z ^ 1 
is given by the counting rules. For this limit, the two 
quarks not in the proton must give their momentum to 
the three that are, and this leads to the factor (1 — z)^. 
Then, barring effects external to the bq residuum, the 
proton should have 3/5 of the residuum's momentum and 
this requires N = 5. There is, however, some pull from 
the struck quark which could increase the probability of 
protons going in the forward direction. This can be ac- 
commodated in the above fragmentation function by re- 
ducing N, and we shall quote results for both iV = 3 and 
N — 5. Lower values of N increase the cited values of /. 

Fig. |l| shows the comparison of the E545 data and 
WA25 data (lO) with the sum of the 6q and 2N con- 
tributions, using N — 5 and / = 2%, or equivalently for 
A^ = 3 and / = 4%. The E545 data is absolutely normal- 
ized, da/dp — {a^^f{i'd)/Ntot)Nbin/Ap where Ap is the 
bin width and a^^f{i'd) is obtained from [p5| and evalu- 
ated at the 50 GeV average E^, in the E545 experiment. 
The BEBC data WA25 is scaled to agree with the E545 
data at the lower momenta. 

The match between the calculation and the data is 
quite good. A 2N contribution alone, with this wave 
function, could not match the data. The 6q contribution, 



though it is overall only 2-4% of the normalization of the 
deuteron state, contributes the major share of the cross 
section for energetic backward protons. The crossover 
momentum is about 300 MeV, and above this momentum 
a larger and larger majority of the protons come from the 
6q cluster, to further elaborate this point, we note that 
at 500 MeV (c = 1) backward proton momentum, the 
Hulthen contribution needs 800% additional contribution 
to be in agreement with the new data. 
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FIG. 1. Protons in backward hemisphere from neutrino in- 
duced deuteron breakup; p is the magnitude of the three mo- 
mentum. The data is Fermilab 15-ft bubble chamber (E545) 
data and BEBC (WA25) data. The 6q curve as shown uses a 
fragmentation function whose transverse part has a power-law 
falloff, (1 + Pt/A^)~^, and whose longitudinal part goes as 



'(1 



and is normalized to give one proton per 6q clus- 



ter breakup. The normalization requires a 2% fraction of 6q 
state. Using instead z'^{l — z)^ leads to a curve whose shape 
is essentially the same as the one shown in the important re- 
gion above 250 MeV/c, and with the same normalization if 
the 6q fraction is 4%. The 2N calculation uses a Hulthen 
wave function pi| , adjusted for use in a light front version of 
a relativistic calculation. 

Backward proton data from neutrino scattering is also 
available from Argonnc EQ], Brookhaven pOl, and again 
Fermilab E545 Il8| for the "quasi-elastic" reaction. 



^ +p + Ps, 



(22) 



where Ps is a label for the "slow protons" or "spectators." 
Fig. shows the backward proton spectra from the 
three sets of "quasi-elastic" data, together with the spec- 
tra from the inelastic data of E545 and the WA25 data. 
There is reasonable consistency within errors among all 
the data sets for the proton spectrum, despite great dif- 
ferences in the incoming neutrino energy. One may think 
that the material struck by the incoming probe goes for- 
ward and that any backwardly emerging hadrons have 
spectra governed only by the distribution of constituents 
in the target. The consistency among the data sets for 
backward protons supports this view. 



To repeat the main point of this section, a calculation 
modeling the deuteron as 2N plus a small amount of 6q is 
able to match the data out to about 720 MeV backward 
hemisphere proton momentum, again within experimen- 
tal uncertainty. 
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FIG. 2. Protons in the backward hemisphere from several 
different neutrino induced deuteron breakup experiments; p 
is the magnitude of the three momentum. The data labeled 
E545 and WA25 are for inclusive deuteron breakup, and the 
data labeled ANL, BNL, and FNAL (the latter actually also 
from E545) are for quasi-elastic deuteron scattering. The ver- 
tical scale is for the E545 inclusive data; the other data are 
scaled to it. 



IV. USE OF OTHER WAVE FUNCTIONS 

The two-nucleon wave function that we have used in 
the above discussion is fairly simple, and one may inquire 
what happens if a more realistic wave function is used. 
There are many wave functions derived from nucleon- 
nucleon potentials that are fit to nucleon-nucleon scatter- 
ing data, and sometimes also to electron-deuteron scat- 
tering data. The assumption is made that only nucleon- 
nucleon, or sometimes only baryon-baryon, degrees of 
freedom are needed. If there are other degrees of free- 
dom present — and that is a crucial question we are try- 
ing to address — they are ignored. But if they exist, their 
effects are present in nature, and they must be included 
in the fitted baryon-baryon potential. And if one cal- 
culates deuteron structure from such a potential, who is 
to be sure if small effects in the (high momentum) tail 
of the wave function are really due to the two nucleons, 
or due to the fitted wave function trying to emulate an- 
other degree of freedom? (We should quote the Paris 
group's remark that "there is no compelling theoretical 
reason to believe the validity of our potential in the region 
r < 0.8 fm. since the short range (SR) part of the inter- 
action is related to exchange of heavier systems and/or to 
effects of subhadronic constituents such as quarks, glu- 
ons, etc." |2q] There apparently would be no physical 



significance in the invention of a nucleon-nucleon poten- 
tial model giving enhanced high momentum components 
by modifing this r < 0.8 fm. region. But, the important 
point is: It is precisely this region that this our new data 
probes.) 
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FIG. 3. Comparison of the Bonn, Hulthen, and Paris wave 
functions. (The 0^ on the vertical axis is a sum of 4>% and (f''h 
for the Paris and Bonn wave functions.) 



Thus, if one speaks of more realistic wave functions in 
the present context, one may object to the phrase "more 
realistic." We shall however take the Paris and Bonn 
wave functions p7| , Pq | as representative of more sophis- 
ticated wave functions and see what happens when we 
use them. A first plot, Fig. ^ shows the Paris, Bonn, 
and Hulthen wave function in momentum space. They 
are essentially the same at low momenta, but at several 
hundred MeV, thanks mainly to the shape and 5.77% 
size of the D-state, the Paris wave function is consider- 
ably larger. Among other wave functions, both the Reid 
wave function [p9| and the relatively new wave function 
of Van Orden, Devine, and Gross |3^ are rather close to 
the Paris wave function. 

Fig. Ha shows what one of the Bonn wave functions, the 
energy independent OBEPQ Q produces for the back- 
ward proton spectrum. Though a more complete and re- 
alistic wave function, the results are not strikingly differ- 
ent than those from the Hulthen wave function. Fig. E|b 
shows a corresponding plot for the Paris wave function. 
All these potential models of the nucleon-nucleon inter- 
action give representations of our new data which are 
well below the data for the high momentum backward 
protons. At a momentum of 500 MeV (c — 1), a con- 
siderable contribution must be added to each to bring 
them close to the data being presented: For the Hulthen, 
as noted above, the additional amount needed is 800%; 
for the Bonn model an increase of 350% is needed (see 
Fig. ^a) and for ths Paris model, an additional 62% is 
needed to bring theory and experiment into agreement. 
The logarithmic scale for da /dp in these figures misleads 
the eye. But, in Fig. 0b, it is clear that the Paris curve 



is below the majority of the data error bars. To make a 
more compelling statement, we give a simple statistical 
comparison of the Paris curves with our data in the mo- 
mentum range, 0.2 < p < 0.5 GeV. The squares of the 
differences between the center of the data points and the 
dashed Paris and solid Paris -|- 6q curves divided by the 
error bar is calculated. The result is a value of x^ per 
data point of 0.8 for the solid curve and and 3 for the 
dashed curve, corresponding to confidence levels of 0.5 
and 0.003, respectively, where the Particle Data Groups, 
graphs are used as a most universal particle physics con- 
vention. 
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FIG. 4. Comparison of the backward proton spectrum in 
v + d ^ i^ + Pb + X to calculations using the Bonn (en- 
ergy independent OBEPQ) wave function and the Paris wave 
function, each using a z^{l — z)^ form in the fragmentation 
function and with a 2% 6q contribution for the Bonn wave 
function case, 1% for the Paris case. As in Fig. hi, similar 
results follow with a 2'^(1 — z)'^ form and a 4% or 2% 6q con- 
tribution for the respective cases. 



V. DISCUSSION AND CONCLUSIONS 

We have seen that a small amount of 6q cluster in a 
deuteron can explain the backward proton data. It is still 
possible that some 2N wave function with an increased 



amount of probability at high momenta could also ex- 
plain the data. Therefore it is of interest to find other 
signatures that could signal the presence of the 2N or 6q 
states. We will mention two possibilities, one if a po- 
larized deuteron target is available, another if there is 
enough data to bin in both x and p, and then conclude. 
If a polarized target is available, then the 2N model 
leads to characteristic variations of the backward pro- 
ton angular distribution. At low momentum, the wave 
function is mostly S-wave and the backward distribu- 
tion is angle independent regardless of the deuteron's 
polarization. At a momentum where the D-state dom- 
inates (about 400 MeV for the Paris wave function, as in 
Fig. H), a polarized deuteron has a non- isotropic spatial 
wave function. Fig. |5| shows the angular distribution of 
backward protons from the D-state for both a longitudi- 
nally polarized deuteron (i.e., Jz = with quantization 
direction along the incoming current) and for the average 
of the two transverse polarizations. (The latter would be 
for either transverse polarization in the analogous elec- 
tromagnetic case.) 
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FIG. 5. Angular distribution of backward protons assum- 
ing a polarized deuteron, the 2N model, and a momentum 
at which the D-state dominates. The solid curve is for the 
average over the two transverse polarizations of the deuteron, 
the dashed curve is for the longitudinal polarization, and the 
dotted hne shows the sum over all polarizations. 

If we can bin data in both x and p, then we can define a 
"two-nucleon test ratio." This is simply the ratio of the 
observed differential cross section for backward proton 
scattering to the cross section for scattering off a neu- 
tron with an appropriate momentum shift. The latter is 
intended to be just what would be expected for the 2N 
model, with the wave function factor removed. Explic- 

itly, 



i?i- 



,{x,y,a,pT) 



anx{x,y,a,pT) 
where the denominator is 
da 



(23) 



CTnX ~ 



dxdy 



{vn -^ jjL X) 



- K{2 - a) [D^iO + 5„(0 + (1 - yfU^iO] ■ (24) 
If the 2N model is correct, then 

i?i = |V'(a,pT)|'. (25) 

Thus, we can test for the 2N model by plotting i?i vs. 
X at fixed a and pt, or at fixed p. If the 2N model is 
right, such a plot would produce just a simple horizontal 
line. The crucial question is how different a result a 6q 
cluster would give. We have elaborated on this question 
in Rcf. p^] for the electromagnetic case. If the 6q cluster 
dominates at some fixed backward proton momentum, 
it gives a curve for i?i vs. x that varies by a factor 
of roughly two from peak to valley. It should be easily 
distinguishable from the 2N expectation. 

In conclusion, we have studied the production of back- 
ward protons in neutrino- and antineutrino-deuteron 
scattering, and compared the existing data to one model. 
The backward proton data from the E545 Fermilab ex- 
periment shown in this paper has not previously been 
published, although it has appeared in talks Q. The 
model we have considered is an incoherent sum of con- 
tributions from 2N and 6q components of the deuteron. 
None of the 2N models that we have looked at has by 
itself a large enough high momentum tail to explain the 
backward proton data above about 300 MeV/c. If we add 
a 6q component, we get straightforwardly a good match 
to the shape of the backward proton spectrum at high 
momentum. If the probability of the 6q cluster is one to 
a few percent in the deuteron, then the 6q contribution 
accounts well for the observed normalization of the data 
at high momentum, while adding negligibly to the 2N 
contribution below about 250 MeV/c p% - 

We consider this a good indication that 6q configura- 
tions exist in the deuteron and can be observed in certain 
circumstances. It is however not ironclad proof since in 
principle it may be possible that some 2N wave function 
with an enhanced high momentum tail could also explain 
all the backward proton data. But one should realize that 
a wave function gotten from a potential that is fit to data, 
including low energy nucleon-nucleon scattering data, is 
matching a Nature that may contain 6q cluster effects 
and must mock them up somehow in the context of its 
own degrees of freedom. This means that a good fit to 
the data with just a 2N wave function is not in its own 
turn ironclad proof against a 6q cluster. Hence we have 
added suggestions of further tests that may eventually 
argue directly against the 2N models. 
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